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By use of external periodic driving sources, we demonstrate the possibility of controlling the 
coherent as well as the decoherent dynamics of a two-level atom placed in a lossy cavity. The con- 
trol of the coherent dynamics is elucidated for the phenomenon of coherent destruction of tunneling 
(CDT), i.e., the coherent dynamics of a driven two-level atom in a quantum superposition state 
can be brought practically to a complete standstill. We study this phenomenon for different initial 
preparations of the two-level atom. We then proceed to investigate the decoherence originating from 
the interaction of the two-level atom with a lossy cavity mode. The loss mechanism is described 
in terms of a microscopic model that couples the cavity mode to a bath of harmonic field modes. 
A suitably tuned external cw-laser field applied to the two-level atom slows down considerably the 
decoherence of the atom. We demonstrate the suppression of decoherence for two opposite initial 
preparations of the atomic state: a quantum superposition state as well as the ground state. These 
findings can be used to the effect of a proficient battling of decoherence in qubit manipulation pro- 
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I. INTRODUCTION 



The idea of controlling the coherent dynamics of a quantum system by an external time-dependent force has found 
wide spread experimental and theoretical interest in many areas of physics (for reviews see fl|-H). It is, e.g., a 
commonly used tool to manipulate trapped atoms in quantum optics Q,^] as well as to control chemical reactions 
by a strong laser field In the context of quantum optics, it has been demonstrated experimentally || that 

a frequency-modulated excitation of a two-level atom by use of a microwave field driving transitions between two 
Rydberg Stark states of potassium significantly modifies the time evolution of the system. In the context of tunneling 
systems it has also been demonstrated that it is in principle possible to completely suppress the coherent tunneling 
of an initially localized wave packet in a double-well potential by an external, suitably designed time-periodic cw- 
perturbation {coherent destruction of tunneling) 0. 

However, real quantum systems are always in contact with their environment. The coherent dynamics is then 
usually destroyed due to the influence of the large number of environmental degrees of freedom. Not only the phase of 
the quantum system is disturbed (decoherence) H but also energy exchange (dissipation) |9|,[l0| takes place between 
the system under consideration and the environment. An example of such a system- bath interaction is the ensemble 
r^j* of electromagnetic field modes in a cavity, each of which is described as a quantum mechanical harmonic oscillator 
j|. Each mode interacts with an atom trapped in the cavity. On the other side, the cavity modes themselves are 
also not isolated from the macroscopic environment; as such they are more realistically described as damped qua ntum 
harmonic oscillators. A topic of fundamental interest is the decay of quantum superpositions of states. In pl[ it is 
shown how quantum optical nonclassical states are highly sensitive to dissipation stemming from a zero-temperature 
heat bath. Experimental works studying decoherence systematically are rare. In [ fl2| the decoherence of mesoscopic 
superpositions of field states in the cavity has been investigated. In a recent work, Wineland and collaborators |l3| ] 
demonstrate that the decoherence rate scales with the square of a quantity that describes the separation between two 
initial states. Moreover, Knight and co-workers ji~4| proposed an experimental scheme to probe the decoherence of a 
macroscopic object. 

In this spirit, the question arises to which extent it is possible to control the dynamics of a quantum system in 
presence of decoherence and, moreover, whether the effect of decoherence can be minimized by an external time- 
dependent force, e.g., by a laser field ||l|,|^,|l5| ^18| , p0| . To achieve this goal, various approaches have been undertaken 
in recent years, (i) It has been shown that the effect of coherent destruction of tunneling (see above) can be used to 
slow down the relaxation of a quantum system to its asymptotic equilibrium J16[. (ii) Moreover, a suitable tailored 
sequence of radio- frequency pulses ("quantum bang-bang" [O] or "parity kicks [|18[) that repeatedly flip the state of 
a two-level atom may be used to suppress decoherence. (iiijThe cavity-induced spontaneous emission of a two-level 
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atom can be manipulated by a strong rf field which couples to the cavity mode p9| . (iv) The manipulation of the 
system- bath interaction by a fast frequency modulation also results in slowing down decoherence and relaxation [^o| . 

The objective of this work is to study the influence of a time-periodic driving field on the dynamics of a two-level 
atom. In the first part of this work (section |lj]), we deal with the objective to "freeze" a coherent dynamics, i.e., 
we shall employ an effect known as coherent destruction of tunneling. Most importantly, we investigate this freezing 
phenomenon from the viewpoint of its dependence on different initial preparations. 



In the second part of this work (section III), we do not elaborate further on the effect of coherent destruction 
tunneling, but instead investigate the control of decoherence of a two-level atom placed in a lossy cavity. Our model 
consists of a two-state system which is coupled to a time-dependent periodic field. The driven two-state system 
interacts furthermore with one mode of the cavity having the frequency fi. This mode is itself damped by the 
coupling to a bath of harmonic oscillators (lossy cavity). It is known ]2l[ | that a Hamiltonian consisting of (1) a 
system part, (2) a harmonic oscillator with frequency f2 that is being coupled to the system, and (3) a bath of 
harmonic oscillators which are coupled to this very single harmonic oscillator can be mapped onto a Hamiltonian 
composed of the system part coupled to a harmonic bath with an effective spectral density. This effective spectral 
density possesses a Lorentzian-shaped peak at f2. The completely isolated atom (no driving, no cavity) evolves in 
time in a coherent way according to the Schrodinger equation. It is this dynamics which we want to preserve and 
protect as far as possible from the decoherent influence of the environment. Our major finding is that a cw-control 
field can indeed be used to (i) reduce decoherence and (ii) to restore to some extent the unperturbed, non-dissipative 
time-evolution. 

II. THE DRIVEN TWO-LEVEL ATOM 
A. Floquet Formalism 

To start we consider a Hamiltonian describing a two-level atom with the ground state |1) and an excited state |2). 
The energy levels are separated by the energy HAq. The atom with the transition dipole moment fi is driven within 
the long wavelength approximation by an external, time-dependent laser field of the form E(t) = £q cos(wLi) with 
frequency lj-^ and amplitude £oj yielding the driven quantum system 

H(t) = ~[A a z + s(t)a x ] . (1) 

Here, the matrices &i, i = x,y, z denote the Pauli spin matrices. The part involving s(t) = s cos(uLi) with s = l^iE^fh 
presents the time-dependent driving which couples to the transition dipole moment /i of the atom. Note that within 
this scaling the amplitude s possesses the dimension of a frequency. The driven time evolution of the populations of 
the energy levels exhibits an oscillatory behaviour. For an initial preparation of the atom in the ground state and 
for a resonant driving, i. e. wl = Aq, and with s not large we can invoke the rotating wave approximation. The 
population of each state then oscillates between and 1 with the Rabi frequency Or = s/2. Because the Hamiltonian 
(|l|) is periodic in time with the period T = 27t/wl, i-e., H(t +T) = H(t), we next apply for the general case away 
from resonance the Floquet formalism |l| . The time-dependent Schrodinger equation may be written as 

{H(t) - iHd/dt}\ip(t)) = 0. (2) 

According to the Floquet theorem, there exist solutions to eq. (||) of the form 

|¥ Q (t)> = exp(-k a f/ft)IM*)> , (3) 
with a — 1,2. The periodic function |<I> Q (i)) are termed the Floquet modes and these obey 

\Z a {t + T)) = \Z a {t)). ( 4 ) 

Here, e a is the so called Floquet characteristic exponent or the quasienergy, which is real-valued and unique up to 
multiples of Htu^. Upon substituting eq. (^) into the Schrodinger equation (g) one obtains the eigenvalue equation for 
the quasienergy e a 

H(t)\$ a (t)) =e„|*„(t)) (5) 

with the Hcrmitian operator 
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H{t) = H{t)-\hd/dt. (6) 

We stress that the Floquet modes 

|$a'(*)) = |*a(<)) exp(inw L t) = |$an(t)) (7) 

with n being an integer number n = 0, ±1, ±2, . . . yield equivalent solutions to eq. (|^) but with the shifted quasienergy 

e a — > e a ; = £ Q + n/iw L = e a „ . (8) 

Therefore, the index a corresponds to a whole class of solutions indexed by a 1 — (a,n). The eigenvalues {e a } can 
thus be mapped into a first Brillouin zone obeying —huj-^/2 < e < Huj-l/2. It is clear that for our choice of the external 
driving force, i.e. s(t) = s cosuj^t the quasienergies are functions of the driving amplitude s and the driving frequency 
wl- For adiabatically vanishing external driving they merge into the eigenvalues of the time- independent part of the 
Hamiltonian (|l|), i.e., 

£ an {s, w L ) — ► T^A /2 + nhuj L , (9) 

where the negative (positive) sign corresponds to a = 1 (a = 2). The Floquet modes, correspondingly, turn into the 
eigenfunctions |a) multiplied by an additional phase factor, i.e., 

\* m (t)) ^ |a) exp(iu; L ni) . (10) 

For a finite driving strength s ^ 0, the determination of the quasienergies e a requires the use of numerical methods. 
The interested reader is referred in this context to the literature (2^j2^]. However, we here state without proof that 
in the high-frequency regime Ao -C max[uL, (scjl) 1 ^ 2 ] the difference between the two quasienergies is given by |l5| 

£2,-1 - £i,i = 7jA Jo(s/wl) , (11) 
where Jq denotes the zeroth-order Bessel function of the first kind. 



B. Freezing the coherent dynamics of a driven two-level system 

Eq. (|ll]) implies a most interesting consequence for a driven two- level system |Tq| : If one chooses the driving 
parameter s and cjl in such a way that the argument of the Bessel function is at a zero of the Bessel function, the 
splitting between the quasienergy vanishes. Possible transitions between the Floquet states are then at most induced 
by the remaining periodic time-dependent parts of the corresponding Floquet modes |$ Q (t)). This effect has been 
discovered in the context of tunneling systems. There, a wave packet being an equally weighted superposition of the 
symmetric and antisymmetric ground state is initially localized at one side of a double-well potential. By applying 
an external suitably tailored periodic field, the wave packet can be stabilized and can be prevented from coherently 
tunneling back and forth between the two wells, i.e., one finds coherent destruction of tunneling (CDT) pLQ]. We 
emphasize here that the crossing of two tunneling related quasienergy levels yields a necessary (but not sufficient) 
criterion for the suppression of coherent tunneling fl5|| . 

The challenge we want to address next is as follows: How does the driven dynamics of a two level atom that is 
being prepared in some arbitrary initial state evolve when the corresponding condition for the parameters obey the 



CDT-condition (11)? The system dynamics can be described by its density operator p(t), which is a 2 x 2-matrix, 
i.e., 

p{t)=I/2+ (12) 

i=x : y ,z 

where the expectation values Oi(t) := Tr{p(t)ai},i = x,y,z are the dynamical quantities of interest. / denotes the 
unit matrix and a x (t) and a y (t) are related to the coherences (the off-diagonal elements) of p(t) while a z (t) is the 
population difference between the two energy eigenstates \a). This implies that the state of the quantum system at 
time t is fully determined by the knowledge of the three expectation values o~i(t). 

To determine the state of the driven two-level system at time t, we consider the Heisenberg equation of motion 
for the density matrix. Using the commutation relations for the &i, we arrive at the equation of motion for the 
expectation values ai{t) in (fL2h, i.e., 
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<Mt) = -A cr.y(i) , 

cj y (t) = A a x (t) - s(t)a z (t) , 

& z (t) = s(t)a y (t) . (13) 

To study the dependence of the effect of coherent destruction of tunneling on the initial preparation we first choose 
as initial state an equally weighted coherent superposition of the two unperturbed energy eigenstates, i.e., 

|*(t = 0)> = -^(|l) + |2)), (14) 



corresponding to a x (t = 0) = l,cr y (t = Q) = a z (t = 0) = 0. We solve the set of coupled differential equations (|1J 
numerically by a standard fourth order Runge-Kutta integration algorithm with adaptive step-size control. In Fig. 
[l]a, the time-dependence of the three expectation values is depicted. The driving parameters are chosen such that the 
condition ( p"l| ) is fulfilled: in doing so we use Cl>l = 50Ao and s = 120. 241. ..Ao- Surprisingly all three expectation 
values <Ji{t) can be brought simultaneously to an almost perfect standstill ! 

Next, we choose the ground state as initial state, i.e., we use \^(t = 0)) = |1). This corresponds to a x (t = 0) = 
o y (t — 0) = 0,<J z (t = 0) = 1. The result is depicted in Fig. [l]b. Applying to the so prepared two-level system a 



laser field obeying the CDT-condition (11), we find that the y-component cr y (t) and the z-component o~ z (t) exhibit 
strong oscillations. This oscillations follow from the numerically evaluated Floquet theory for the driven two-level 
system, and are not described by the Rabi-oscillations as predicted from a rotating wave approximation; this latter 
approximation is strongly violated for our chosen set of driving parameters. In contrast, o~ x (t) can be stabilized 
around the initial value of zero. This finding is in accordance with the CDT-phcnomcnon: it reflects the fact that 
the corresponding two equally weighted ("left" and "right") localized parts of the ground state wave function of a 
double-well potential, as represented within a localized representation, can each be stabilized too. 

This CDT-effect opens the doorway to manipulate the influence of an environment on a quantum system. It 
is known ]lq| that the coherent destruction of tunneling survives to some extent in presence of a coupling to the 
environment. Certainly, the system will relax in presence of an environment; however, as it is demonstrated in p6| , 
the relaxation process can considerably be slowed down in the presence of a CDT-ficld. 

In view of using differing initial preparations, the following remark should be made. From the viewpoint of stabilizing 
the state of a qubit (characterized by a quantum mechanical two-level system) in a quantum information processor 
|p4j , it is of foremost interest to stabilize the coherent superposition of two states of the qubit. Thus, our first 
choice (|l4|) is of relevance in the context of the possibility for quantum computing. Moreover, fundamental questions 
concerning the decoherence of superposition states arise for the physics that occurs when one crosses the interface 
between the classical and quantum world, and vice versa PJl3]| . 



III. CONTROL OF DECOHERENCE FOR A TWO-LEVEL ATOM 



In this section we shall study the influence of an applied cw-control field for reducing decoherence of a two-level 
atom placed in a lossy cavity. 



A. Driven two-level atom in a lossy cavity 



To start we consider a two-level atom in a dissipative environment, e.g., a lossy cavity wherein the leakage of photons 
damps the radiation field. Additionally, the atom may be manipulated by a time-dependent external field like a laser 
beam. In our model, the driven two- level atom is represented by the Hamiltonian ([!]). It is coupled to one mode 
of the cavity which is described by one harmonic oscillator with frequency f2, characterized by the annihilation and 
creation operators B and W which fulfill the usual commutation relations for bosonic field operators. The coupling 
constant is denoted by g and has the dimension of a frequency. This cavity mode is damped by a bilinear coupling 
to a bath of harmonic oscillators of frequencies u>i. They are similarly described by bosonic annihilation and creation 
operators bi and b\. The coupling constants of the cavity mode to the harmonic bath are given by Ki and have the 
dimension of a frequency. The total system-bath Hamiltonian is therefore written as 

H{£) = --[A a z + s(t)a x ] 

+ hVL{B^B +\) + + B)a x 
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J2 fuji{b\bi + h + h{& + B)J2 Ki(f>l + U) • (15) 



i=l i=l 

The influence of the bath on the two-level atom plus cavity mode is fully characterized by the spectral density 

N 



J(uj) = 2nJ2 K i S ( UJ - UJ i)- (16) 
i=i 

We let the number of bath modes going to infinity (N — » oo) and choose an Ohmic spectral density for the bath 
oscillators with an exponential cut-off at some large frequency uj c 3> Ao,cjl,^, i.e., 

2r 

J{u) = —ujexp(-uj/uj c ) , (17) 

Here, we have introduced the damping constant T which is related to the quality factor of the cavity. Since the cavity 
mode as well as the bath oscillators are described by harmonic oscillators, we follow the approach in |pl[| and map 
the Hamiltonian ( |l5| ) onto a Hamiltonian where the central system, i.e., the two- level atom, is now bilinearly coupled 
to a bath of mutually non-interacting harmonic oscillators with an effective spectral density J b r(uj). Upon letting the 
cut-off frequency going to infinity, i.e., u> c — > oo, this effective spectral density emerges as 

_ ier g^n 2 

JcffM-— (f7 2_ cj2)2+4cj2r2 - (18) 

For small frequencies lo, it increases linearly like in the original Ohmic spectral density J(u>). However, it has a 
Lorentzian shaped peak at ui = Q with a line width r < Q. 

In the following section, we make extensive use of the bath autocorrelation function A4(t) — Ai'(t) + iA4"(t), which 
is obtained in terms of the effective spectral density J e ff(w), i. e. 



1 f°° 

M(t) = — dujj c s(uj) 
t Jo 



coth ^ 2fc J 1 ^ < " s '"' 1 ' s ' nl ~' ! 



(19) 



All our further considerations treat the case when the bath is at zero temperature, i.e. T = 0. In this limit, and for 
our choice of the effective spectral density (fl8|), we obtain for the real and imaginary part, respectively, the analytical 
results 



M '(t) = ^-g 2 tf \~ - 1 e- rt co S (Vn 2 -T 2 t) ~ ( dy 



( y 2 + n 2 ) 2 - vr 2J ' 



M"(i) = -4/-=_= e" ri sin( Vf^Ht) . (20) 

The quantity of interest is the reduced density matrix for the two-level system which we denote - just as in the 
undamped case - by p{t). It follows by tracing over the bath degrees of freedom in the full density operator W(t) 
which corresponds to the system-plus-bath Hamiltonian (|l5|), i.e., p(t) — treW(t). Like in the deterministic case in 
cq. (p^), p{t) is fully characterized by the expectation values <7i(t), i = x,y, z. We shall determine their corresponding 
equations of motions next. 



B. Bloch - Redfield Formalism 

To deal with quantum dissipative systems, several techniques have been developed |l]|^,[l(|^2| . A very efficient 
numerical algorithm for general quantum system with a discrete eigenvalue spectrum has been developed by Makri 
and Makarov within the real-time path-integral formalism [p5[ . It has also been applied to spatially continuous 
tunneling systems in presence of driving (26). Moreover, the real-time path-integral formalism has extensively been 
used to describe a moderate-to-strong (!) two- level system- bath interaction ]j],p|, pl|j27| . Recently, the former scheme 
has been generalized to describe multi-level, driven vibrational and tunneling dynamics in |2§f| . At weak system-bath 
coupling the Nakajima-Zwanzig projector operator theory p9| provides a powerful tool to describe the corresponding 
reduced density matrix dynamics. 

For our quantum optical problem at hand, the suitable method of choice in presence of a physically realistic weak 
system-bath coupling is the projection operator technique: it yields in Born approximation the generalized master 
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equation. It can be simplified further without loss of accuracy in leading order in the (weak) coupling strength g 
by invoking the Markovian approximation |j32| . For a strong harmonic driving this objective was formally (only) 
developed a long time ago by Argyres and Kelley |30|. Following the reasoning in |}Cj (see in this context also 
we recently have derived for this case of a driven spin-boson problem with an arbitrary control field the explicit set of 
coupled, Bloch-Redfield type equations Q 

0x(t) = -A a y (t) , 

& v (t) = A a x {t) - s{t)a z {t) - ri (t)oy(t) - T 2 {t)a x (t) - A v {t) , 

<j,(t) = s{t)a y {t) - Ti(t)o z {t) - T 3 {t)a x {t) - A z (t) . (21) 

The time-dependent rates Ti(t) — J*dt'Ai'(t — t')bi(t,t') , together with the inhomogeneities A y (t) = KeF(t), 
A z (t) = 1xaF(t), with F(t) = (1/2) J* dt'M"(t - t')[u 2 {t, t') ~ v 2 (t, t')] determine the dissipative action of the thermal 
bath on the two-level atom. The functions M' and M." denote the real part and imaginary part, respectively, of 
the correlation function M given in eqns. @. The quantities u(t,t') = (l\U(t,t')\l) + (2\U(t,t')\l) and v(t,f) = 
(l\U(t,t')\2) + (2\U(t,t')\2) are sums of matrix elements of the time evolution operator U(t,t') of the isolated (i.e., 
g = 0) driven two-level system. The functions bi read &i = Reuv* , b 2 = —(1/2) Im(u 2 — v 2 ), and 6 3 = (1/2) Re(u 2 — v 2 ). 
Note that this set of equations is valid in theparameter regime g <C Aq/2. One can demonstrate that for the undriven 
case, i.e., s = 0, the analytic solution of eq. ( |2l| ) in first order in g reproduces the analytical path integral weak-damping 
results in Refs. IQ. 



C. Controlling the decoherence of a quantum superposition of states 



The idea of controlling the decohering influence of the environment on a quantum system by an external time- 
dependent field is demonstrated for the case of the two-level atom which is initially prepared in an equally weighted 
superposition of the two energy eigenstates given by o~ x (t = 0) = l,o~ y (t = 0) = a z (t = 0) = 0. In doing so, we 
consider four different situations: (1) first, we look at the isolated two-level atom dynamics without driving and 
without coupling the atom to the lossy cavity mode. This case corresponds to setting s — and g = 0. Case (2) is 
devoted to the driven two-level dynamics. We switch on a coherent driving cw- field but keep the system isolated from 
the bath, i.e. s / and g — 0. In case (3) we investigate how the undriven system dynamics relaxes in presence of a 
dissipative coupling to the bath. We therefore set s = and g ^ 0, T ^ 0. Finally, we demonstrate with case (4) how 
this decoherent dynamics can be manipulated with the help of an externally applied time-dependent control field and 
set s ^ 0,g ^ and T ^ 0. 

In order to preserve the coherent evolution of the two-level atom and to protect it as far as possible from the deco- 
herent influence of the environment, we choose the following control scheme: Guided by the physics of a rotating wave 
approximation for the driven system that most closely retains the unperturbed dynamics of an initial superposition 
state (Q) we choose the following parameters: The frequency and the amplitude of the driving field are taken to be in 
resonance with the level spacing of the two-level system, i.e., wl = Ao and s = Aq which corresponds to a moderately 
strong driving strength. This choice implies for the ratio of the corresponding Rabi-frequency and driving strength the 
value 0.5. This indicates that the rotating wave approximation should be used already with caution. Note that under 
the CDT-condition in (11), the field strength would assume an even larger value of s — 2.4048Aq. For the strength 
of the coupling between the two- level atom and the cavity mode we assume g — 0.05Ao. This value is consistent with 
the range of validity of the Bloch-Redfield formalism in Born approximation (see above). The dissipative system-bath 
mechanism is specified as follows: the frequency of the cavity mode is chosen to be in resonance as well, i.e, 57 = Ao. 
By doing so, we in essence maximize the influence of the bath. For the line width of the cavity mode we set T = 0.1 Ao. 
This rather large value mimics (on purpose) an extreme situation because the line width in most realistic situations is 
in general much smaller: Nevertheless, such smaller values would intensify our appealing finding of a driving-induced, 
enhanced recovery of coherence even more. Moreover, the temperature is always set to T — 0. 

Our novel results are depicted in the Fig. 2 (a)-(c) and the Fig. 3 (a)-(c). Fig. ||a depicts the time-evolution of the 
x-component o~ x (t). The isolated two-level dynamics (dashed line) shows coherent oscillations between -1 and 1 at the 
frequency of the level spacing Ao. On top of this line one finds (barely visible dotted line) the results for the driven two- 
level dynamics. This good agreement follows also from the corresponding rotating wave approximation, yielding for this 
preparation just the undriven result. The decoherence in presence of a finite bath coupling (g = 0.05Ao,r = O.lAo), 
see the dashed-dotted line, yields an oscillatory decay towards equilibrium o~ x (t — > oo) = 0, whose envelope is made 
visible by the connecting solid line. Next we switch on the cw-laser control field. As a main result we find that 
the decoherence becomes considerably slowed down - following closely the isolated driven dynamics. This enhanced 
recovery of coherence for the dissipative driven dynamics is made visible to the eye by the connecting weakly decaying 
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and oscillating envelope. This surprising result is rooted in the following facts: The dissipative, non-driven dynamics 
experiences a most effective dissipation. This is due to the resonant coupling at Q = Ao of the two-level atom to bath 
with the effective spectral density in ( |l8| ) which peaks at u = Q. In contrast, the strong driving now dresses this level 
spacing, and moves it out of resonance with the lossy cavity mode. This results in a considerable slow down of driven 
decoherence for <r x (t). 

The decoherent dynamics of the y-component a~ y (t) is qualitatively similar to a x (t). It is depicted in Fig. ||b for 
the same choice of parameters. 

The population difference cr z (t) is shown in Fig. ||c. For the isolated two-level dynamics, o~ z (t) remains constant 
at zero (dashed line) since the system is in an equally weighted superposition of two eigenstates, yielding an obvious 
zero population difference. In presence of the cw-laser control field, the driven dynamics (dotted line) yields a finite 
oscillation of population difference. This deviation from zero also reflects the deviation from the corresponding 
rotating wave solution (being identically zero for this preparation). Nevertheless, this driven dynamics still exhibits 
an approximate periodicity that closely coincides with the Rabi- value Qr = Aq/2. 

The undriven, dissipative relaxation to equilibrium (dashed-dotted line) proceeds with temperature T = almost 
completely towards the ground state with corresponding maximal population difference o~ z (t — > oo) ~ 1. Due to 
the coupling to the cavity mode performing zero point oscillations, the value of 1 is not fully reached. The driven, 
dissipative relaxation (solid line) to the time periodic asymptotic state exhibits oscillations around zero - following 
initially (up to Aot w 50) closely the driven coherent dynamics. In virtue of Floquet theory for the long-time limit 
of the time-periodic generalized Bloch-Redfield equations in (11), this asymptotic periodicity matches in the lon-time 
limit the frequency of driving, i.e., u)l = Aq (not depicted). 



D. Controlling the decoherence from the atom ground state 



To answer the question whether the proposed control scheme works as well in the opposite limit of an initial state 
which is an eigenstate we next choose the ground state as the initial preparation, i.e., we use a x (t — 0) — a y (t = 0) = 
0, a z (t = 0) = 1. The remaining parameters are taken to be the same as in the previous subsection III C. 

Fig. ||a shows the decoherent dynamics for cr x (t). Since the chosen initial state is an eigenstate of the isolated two- 
level system no dynamics is exhibited (note the filled squares on the line at zero in the figure) . This situation remains 
unaltered in presence of a dissipative coupling of the quantum system, as indicated by the asterisks on the line at zero. 
At zero temperature the system at weak dissipation remains essentially in its ground state. Upon switching on the 
driving with no coupling to the lossy cavity present, the driven two-level dynamics exhibits a Rabi-like quasiperiodic, 
oscillatory behaviour (dotted line) . This nonperiodic behavior is rooted in the deviation of the full Floquet dynamics 
from a rotating wave prediction. With our strong driving strength we a priori cannot expect good agreement with 
the corresponding rotating wave approximation. The coupling to the lossy cavity mode damps this quasiperiodic 
behaviour, following for short times the driven isolated dynamics (see solid line), before settling down to asymptotic, 
long-time oscillations at the frequency of driving wl = Ao with a finite, but strongly reduced amplitude (not depicted). 

The decoherent dynamics of the y-component cr y (t) is again qualitatively similar to that of o~ x (t). It is presented in 
Fig. ||b for the same set of coupling and driving parameters. 

Finally, the time evolution of the population difference o~ z (t) is depicted with Fig. ||c. Clearly, the isolated dynamics 
from a prepared initial ground state remains constant at o~ z (t) = 1 (filled squares). The driven dynamics of the two- 
level system exhibits strong non-detuned Rabi oscillations at frequency flu = Aq/2 between -1 and 1 (dotted line). In 
this case the rotatating wave prediction (not depicted) actually yields surprisingly good qualitative agreement with 
the exact dynamics. 

The case of no driving (s = 0) but with a coupling to the bath (g — 0.05Ao,r = O.lAo) shows again a trivial 
dynamics. It relaxes in this case of weak dissipation with a small relaxation rate towards a slightly reduced constant 
value close to 1 (indicated by the asterisks). 

The case with resonant driving (s = Ao,^l = Ao) switched on and simultaneous coupling to the lossy cavity 
mode (with g — 0.05A ,r = 0.1A ) exhibits damped Rabi-oscillations (solid line); it eventually settles down in the 
asymptotic long-time limit to periodic asymptotic oscillations at twice the Rabi frequency and amplitude smaller than 
1 (not depicted). 



IV. CONCLUSIONS 



In this work we have investigated the possibility to control the time-evolution of a two-level atom by time-dependent 
external, periodic control forces. We have demonstrated that the coherent dynamics of the system can be brought 
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to an almost perfect standstill by choosing the ratio of driving amplitude s and driving frequency cjl at a zero of 
the Besscl function Jo(s/wl) (coherent destruction of tunneling). For an initially prepared quantum superposition 
of states all three components oi,i = x,y,z and therefore the entire density matrix p can be locked simultaneously. 
For the initially prepared ground-state, the x-component a x can be stabilized; the other two components a y and a z , 
however, depict strong (non-Rabi) oscillations. 

In presence of decoherence in a lossy cavity we illustrate that the atomic states can be dressed by a time-dependent 
force which moves the atom and the cavity mode out of resonance. As a consequence, decoherence becomes strongly 
suppressed. We have illustrated this effect for two different initial preparations of the atom: (i) for a quantum 
superposition of states we show that the decoherence can be suppressed efficiently, (ii) The second preparation uses 
the ground-state wave function of the isolated system. In that case the decoherence may also be slowed down, but 
the decohering dynamics never approaches again the initial state. 

These findings put the idea across that the method can be used to bring back the state of the atom close to its 
initial preparation. For the case (i) of a superposition state as initial state the decoherent dynamics of the x- and y- 
componcnt a x , cr y are similar to the undriven dynamics of the isolated two-level system (qubit) . Even the ^-component 
o~ z of the driven dissipative dynamics matches at distinct instants of time the undriven non-dissipative dynamics. For 
the second case (ii) of the ground state as initial state this idea, however, seems to fail for the z-component a z . 

To summarize, our proposed scheme for controlling the coherent and decoherent dynamics of a two-level atom works 
very well for initially prepared quantum superpositions of states. This presents good news for the manipulations of 
quantum bits (two-level systems) being in a superposition of states. It is this very feature which makes quantum 
computation interesting and superior to classical computation. 
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FIGURES 
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FIG. 1. Fig. hp,: Time- dependence of the expectation values aj(t),j = x,y,z, determining the dynamics of the full den- 
sity matrix of the two-level atom according to ( |l2[|l3| ). The driving parameters are chosen such that the condition 
for coherent destruction of tunneling is fulfilled, i.e., the driving amplitude is set to s = 120.241Ao and the driving fre- 
quency u>l = 50Ao. The two-level system is prepared in an equally weighted superposition of the two energy eigenstate, i.e., 
cr x (t = 0) = 1, Oy{t = 0) = <Tz(t = 0) = 0. Fig. [ij): The same as in Fig. |l|a but for an initial state preparation being the ground 
state of the two-level system, i.e., a x (t = 0) = a y (t = 0) = 0, a z (t = 0) = 1. 
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FIG. 2. Time-dependence of the expectation values a x (t) (Fig. |2ja), cr y {i) (Fig. ^b) and a z (t) (Fig. ^|c) for the two- level atom 
initially prepared in an equally weighted superposition of the energy eigenstates, i.e. a x (t = 0) = 1, a y (t = 0) = <r z (t = 0) = 0. 
Shown are four cases: (1) no driving (s = 0), zero system-cavity mode coupling (g — 0) (dashed line), (2) with reso- 
nant driving (s = Ao,wl = Ao), but zero coupling (g — 0) (dotted line), (3) zero driving (s = 0), but with finite cou- 
pling (g = 0.05Ao,F = O.lAo) (dashed-dotted line) and (4) with resonant driving (s = Ao,u>l ~ Ao) and with coupling 
(g = 0.05Ao, r = O.lAo) (full line). The temperature is chosen to be T = and the cavity-mode frequency is set to Q — Ao. 
As a guide for the eye, we mark the envelope of the decaying oscillations by solid lines in Fig. |^a and Fig. gb. 
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(a) a s (t = 0) = (b) a y {t = 0) = 
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FIG. 3. Time-dependence of the expectation values a m (t) (Fig. c y (t) (Fig. and <r z (i) (Fig. ^c) for the two- level 
atom initially prepared in the ground-state, i.e., a x (t = 0) = a y (t = 0) = 0,a z (t = 0) = 1. Similar to Fig. || four 
cases are depicted: (1) no driving (s = 0), isolated two-level system (g = 0) (dashed line with filled squares ■), (2) 
with resonant driving (s = Ao,u>l = Ao), no dissipation (g = 0) (dotted line), (3) no driving (s = 0), with dissipation 
(g = 0.05Ao,r = O.lAo) (dashed-dotted line with asterisks *) and (4) with resonant driving (s = Ao,u>l = Ao) and with 
dissipation (g = 0.05Ao,r = O.lAo) (full line). The temperature is chosen to be T — and and the cavity-mode frequency is 
set to 11 = Ao. 
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